Radiative Leptonic Decays of the charged B and D Mesons Including Long-Distance 

Contribution 
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In this work we study the radiative leptonic decays ol B~~ , D~ and — > jW, including both 
the short-distance and long-distance contributions. The short-distance contribution is calculated by 
using the relativistic quark model, where the bound state wave function we used is that obtained in 
the relativistic potential model. The long-distance contribution is estimated by using vector meson 
dominance model. 
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I INTRODUCTION 

The mechanism of heavy meson decays is one of the 
most interesting and challenging fields in particle physics, 
it involves both strong and weak interactions. Nowadays 
strong interaction in the non-perturbative region is still 
an unsolved problem. Compared with hadronic decays, 
leptonic decay is simpler. Strong interaction only oc- 
curs within the initial particle. Pure-leptonic decay of 
heavy meson can be used to determine the decay con- 
stant, which describes the possibility-amplitude for the 
quark-untiquark emerging at the same point. The pure- 
leptonic decay is helicity suppressed. The decay branch- 
ing ratio of a pseudoscalar meson P with quark content 
Qq within the standard model is 



B(P^-W) = 



g f\ v Qi 
8tt 



- Tp f pm \m P (l-^\)\ (1) 



where Gf is the Fermi coupling constant, Vq q the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix element, 
Tp the life time of the meson P, and mp and mi the 
masses of the meson P and lepton I, respectively. The 
decay rate is proportional to the lepton mass squared mf 
is the consequence of the helicity suppression. However, 
the presence of one photon in the final state can compen- 
sate the helicity suppression. As a result, the radiative 
leptonic decay can be as large as, or even larger than 
the pure-leptonic decay mode. It thus opens a window 
for detecting the dynamics of strong interaction in the 
heavy meson or studying the effect of strong interaction 
in the decay. 

The radiative leptonic decay rates of the charged B and 
D mesons have been studied with various methods in the 
literature. In Ref. [J, B and D s are calculated 

in a non-relativistic quark model, the branching ratios of 
the order of 10 -4 for D s — > Ipj and 10 -6 for B — > 
are found. In Ref. 0] with perturbative QCD approach, 



it is found that the branching ratio of 



v^j is of 



the order of 10~ 3 and D + — > e + v^f of the order of 10~ 4 , 
while the branching ratio of B + — > e + vj is at the order 
of 10~ 6 . On the other hand, a smaller branching ratio 
is obtained for Dr s \ — > within the light front quark 



model [3|. Smaller result for £)( s ) — > W"/ is also obtained 
in Ref. |4| within the non-relativistic constituent quark 
model, which gives that the branching ratio of D~ — » W^y 
is of the order of 10~ 6 and D~ — > of the order of 
10~ 5 . The problem of factorization in QCD for B —> Ivy 
is studied in Ref. 

In this work, we study the radiative leptonic decays of 
the charged B, D and D s mesons to Zz>7 including both 
the short and long-distance contributions. The short- 
distance contribution is considered at tree level. The 
wave function of the heavy meson used here is obtained 
in the relativistic potential model previously 0. The 
long-distance contribution is estimated by u sing the idea 
of the vector meson dominance (VMD) |7Hll| followed 
by the transition of the vector meson to a photon. We 
find that the long-distance contribution can enhance the 
decay rates seriously. 

The remaining part of this paper is organized as fol- 
lows. In Sec. II, we present the short-distance amplitude. 
In Sec. Ill, the long-distance contribution is considered. 
The numerical results and discussion are given in Sec. IV. 
Sec.V is a brief summary. 



II The Short-Distance Contribution 



We use P to denote the pseudoscalar meson which is 
composed of a heavy anti-quark Q and a light quark q, 
such as B and D mesons. There are four Feynman di- 
agrams contributing to the radiative decays P~ 
at tree level, which are shown in Fig. [TJ However the 
contribution of Fig. Q] (d) is suppressed by a factor of 
1 /M^ , it can be neglected for simplicity. The effective 
Hamiltonians corresponding to the other three diagrams 
in Fig. [T] can be written as: 



2 





FIG. 1: Feynman diagrams of short-distance contribution at 
tree level. 



T~La 



-ieG F V Qq 
-ieGpVQq 



\ ~ ?q + m Q 
2 (pyPq) 



P£q(lP L ^), 



m a 



n b = ~~-^ Q q QPt * £q (IP L ^) , 

2(p Y P q ) 



u c = 



V2 
-ieGFVQq 
V2 



QgQP^q U 



" 7 1 ri 



mi 



2(p 7 -P;) 



-Pl^v 



(2) 



where P? is defined as 7 M (1 — 75), and Vq q represents for 
the CKM matrix elements. Qq and are the electric 
charges of the quarks Q and (7, respectively. A is the 
electro-magnetic field. 

1-l a , T-Lb and He can be divided into two terms for con- 
venience, according to the numerator of the fermion prop- 
agator. For example, Ti. a can be written as 

H a = Z!fggg^gg (M al + 2M a2 f (IPl^) , (3) 



2^2 



where 



M» 1 = Q4-!l-<f(l-'y B )q. 
PYPQ 

= - ie ^A a p^v aa +p» (A-v a ) - A" ( Pr v a ) 

1 _ mo — #>„ 
M» a2 = -QA ^7" (1 - 75) q = -A a t^, 

1 P-y'PQ 



(4) 



with 



< = Q 7^(1-75)9, 

P-y-PQ 



(5) 



The amplitude of the radiative leptonic decay can be ob- 
tained by inserting the operator of the effective Hamil- 
tonian between the initial and final particle states. For 



example, the contribution of Fig. [T] (a) is 

A a =< jpi\n a \p> 

(uiPl^Vp) 



-ieQ Q G F V Qq 

~^T2 VU " L '" 



X (< 7 I M a l I P > +2 < 7 I Ma2 I P >) 
(uiPh^Vp) 



-ieQQGFVQq 



2V2 

x (_z e Q ^e; a p 7/ j < I iw I P > +p»e*- < | v a \ P > 
-e»*p r <0\v a \P> ~2e; a <0\t^\P>). 

(6) 



The matrix elements < 



P > and < I tT 



P > only depend on the momenta p p and p 7 . According 
to their Lorentz structure, they can be decomposed as a 
linear combination of two terms of pp and p 7 



Va = < I 



f 



pypq 



-7 M (i-7s)g|P> 



=< 



m P a a m P „ 
AalPp + B al ^p^ 

pp-p-t \ Pp-Pj 

7 M (l-75)9 I P> 



PQ-Pl 



mp 
Pp-Pi 



A a 2PpPp 



pp-p-f 



{Ba2P P P^ 



Ca2P P P«) 



(pp-P-rY 



D a2P "p>i + E a2 m 2 P g afl 



+F a2 ^^€ a ^p Pp p ia 

PP-P-y 



(7) 



The coefficients A al , B al , A a2 , B a2 , C a2 , D a2 , E a2 and 
F o2 are all dimensionless constants. The terms of B a ±, 
C a2 and D a2 do not contribute to the decay amplitude 
A a when substituting the above decomposition into eq|6l 
Therefore these terms can be dropped. The coefficients 
can be obtained by the treatment in the following. Mul- 
tiplying V£ with Pjfj,, we can obtain A a i as 



Aai = — < I Q—^—f 

m P PrPQ 1 



(1-75)9 I P>. (8) 



Similarly, multiplying T" M with p 7Q , we have 

B a .2 + E a2 = 0, 

1 - (9) 

Aa2 = — <Q\Qf p (l- l5 )q\P> . 

Tfl p 

Multiplying with p Ju pp a and g QM , and using B a 2 = 
—E a 2, one can get 



E„ 



B n 



1 



2mp 



<0\Q 



(Pp ■ P 7 ) 



(pp ■pq) 



(1 - 75) 

Pq-P-i 



(10) 



q\P> 
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Finally F a2 can obtained by multiplying T" M with 



in the following 

-ieGpVQg 



Aa-^h+c 



F a2 — 



2m 3 D 



{iVe a ^p lP e; aPPa 



<0| q!^q^ (1 _ 7s)g | P> . (ii) 

PQ ■ Pj 



2V2 

-A [{pp-p^e^* -{ PP -e)p^\ 
p P -e pi-e 



The amplitude At can be treated in the same way with 
some coefficients defined as follows 



1 



(Q q - Qq) 



Pp-Pj PrP-r 



A a2 m D p P 



(14) 



<0|Q l»{l-l 5 )q\P> 

P 7 -Pq 

™2 



mp 
PP - Pi 



Pp ' P~t 



-Pi 



This equation clearly shows that the sum of the contri- 
butions of all the diagrams in Fig. [1] is gauge invariant. 
In eq. ffT4|) the factors V and A are 



<0|Q^ 7 ^(l-7 5 )g|P>= 

Pl'Pq 



mp 
Pp-P~t 



ib2P P Pp 



Pp-Pj 



(B b2P ppi; + C b2 p P p«) 



111 , 



(pp-PiY 



D b2 p"p>t + E b2 mpg° 



V 



A 



(12) 



Q Q m P A al _ Q q m P A bl m P A a2 

PP-Pj PP'Pi PrP-y 

2(Q q F b2 -Q Q F bl )mp 

(pp-P-tf 

Q Q m P A al Q q m P A bl m P A a2 

Pp-P~t PP-P-, 
2 (Q q E b2 - Q Q E b i) mp 



(15) 



Pl'P-y 



Using the matrix element A a2 mpp P —< 
Q7 M (1 — 75) q I P >, the amplitude A c can be treated 
simlarly. Finally, the total amplitude can be expressed 
as 



(PP'Pj) 

Next we shall calculate the coefficients A a \, A a2 , E a2 and 
F a2 - 

The pseudoscalar meson state can be written in terms 
of the quark- antiquark creation and annihilation opera- 
tors 



2V2 L PP-P~f 

^m P A a i 2Q Q E a2 mp 



PP ■ P 7 {pp ■ P 7 Y 



" 2 ( Q<3 ) A a2 m P p p ^ , 



[(pp. Pl )e^ -(p P -e)p»\ 



|p(p = o)>=-L E 



d 3 kV (\k\) 



A b 



Pp-Pi 
-ieGpVQg 
2V2 



V2 



6 l +(fc,tK + (-fc,|) - 6 l +(fc,|)C(-fc,t) I >, 

(16) 



(uiP^Vv, 



Q q m D A b i 2Q q E b2 mp 
Pp-Pi (pp-p-yf 



+2 I Q q PP e ) A a2 m P p p 



Q q m P A bl . a 

— ie ^ P-tpe^ppa 

Pp-P-i 7 

[{PP-P-,)S»*-{PP-E)P»} 



PP-Pj 



A c = 



-ieGpVQg J mpA^ fcQ .„„ 



2V2 



where i is the color index, the factor l/\/3 the normal- 
ization factor for color indices, and k the 3-momcntum 
of the quarks in the rest frame of the heavy meson. The 
wave function ^od^l) has been calculated in the relativis- 
tic potential model previously, the numerical solution of 
the wave function can be fitted in the exponential form 

SE3 



pi-pj 



ie"^"p^E~p Pa 



m P A a2 



Pi -P~t 

Pl'S 



[{pp- Pl )e»* -(pp-e)pfl 



*o \ \k\ ) = 47rVA 3 w pi 



»-A|fc| 



(17) 



PfP-r 



A a2 m P p p 



(13) 



The above equations show that the contribution of each 
diagram in Fig. Q] is not gauge invariant separately, but 
the sum of them is indeed gauge invariant, which is given 



The numerical solutions of the parameter A for D, D s 
and B mesons are quoted from Ref. Q recently 

X D = 3.4GeV _1 , \ Ds = 3.2GcV~\ 
A B = 2.8GeV _1 . 

With the meson state given in eq. (|16l) . the matrix el- 
ement < j QTq j P > can be calculated straightfor- 



4 



rvww 
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FIG. 2: The LD contribution. 



wardly, the result is 



< | QTq I P >= x j ~- / dkcMk^o (k) 



3 1 



(18) 



where M = w g (fc, t)«Q(— fc, |) — u q (k, 1)vq(— k, f) is the 
Dirac spinner for the pseudoscalar meson, it can be ob- 
tained in the Dirac-Representation as 

„_ -*fc +m «) (j+^fe+^y (ig) 

2^m° Q m°(p Q +m Q ) (p° q + m q ) 



III The Long-Distance Contribution 

In this section, we estimate the contributions of long- 
distance physics. According to the spirt of vector meson 
dominance model 04H|, we consider the resonance pro- 
cess P~ — > IvV — > Ivy, where the intermediate vector 
resonance V can be p, lo and <p. The contribution comes 
from the semileptonic intermediate IvV , followed by the 
vector resonance turning to an on-shell photon V — > 7, 
which is shown in Fig. [2] The amplitude of the long- 
distance contribution can be written as 



A 



eG F VQ q Q q 

LD = J1 H - UlT I 1 - 75) VP 



The vector's decay constant /y can be derived from the 
decay rate of V — >• e + e~. After a short calculation, the 
decay constants can be related to the vector meson's lep- 
tonic decay widths 



j-2 P 

J p — 



f = 



fd> — 



4 ™ 2 (Qu - Qd 

3m 3 u 2 
47ra2 (Qu + Qd 
3ml 1 „ 



2 p^ete - 1 



2 oj-^e"re _ 7 



(23) 



Ana 2 Ql 



where a is the electromagnetic fine structure constant, 
Qu, Qd and Q s are the charges of the quarks it, d and s, 
respectively. 

The hadronic matrix element < V | J M | P > in 
eq. (|2"0|) can be decomposed according to its Lorentz 



structure as 



null 



< V I Jv—Ap. I P >= : 

mp + my 

+ i {e V u (m P + m v ) A x (q 2 ) 

£v ■ q 



mp + my 
e v ■ q 



(pp+pv) u A 2 (q 2 ) 



2m P q tl A 3 (q 2 ) [• + i £V 2 q 2m v q u A (q 2 ) 



(24) 



where q = pp — pv, and V , A Q , A%, A 2 and A 3 are form 
factors. 

To obtain the amplitude gauge invariant, we take the 
trick used in Ref.[15| in treating the long-distance con- 
tribution in the precess b — > sy via the resonance J/^>. 
With the Lorentz decomposition of the hadronic matrix 
element in ea. (l24|) . the product of the two matrix ele- 
ments in eq. ([2H|) can be calculated to be 



E- 

y 



£"* < I qj a q I V > 



Py - m v + im v T v 



x <V \Q^(l- l5 )q\P> 



(20) 



For the decays of D and B mesons, V represents for p 
and u mesons, while for D s meson decay, V represents 
for (f> meson. 0y is the relative phase between the long 
and short-distance contributions. The matrix element 
< I qy^q I V > is used to define the decay constant of 
the vector meson 



<0\q^q\V >= C qV fve v \ 
where the factor C g y's are 



(21) 



£ * a < I qi a q |V><7| Jy_ A I P >= 



eQqCqvfv 



-iAi (q 2 



-i (pp-e^ 



mp + my 



e a ^ P PaPvpe; a V (q 2 



^ (m P + m v ) -Py 



pp-pv 



(m P + m v )py 



pp-pv 



Ax (q 2 



(Pi 



Pv) 



mp + my 



■A 2 (q 2 



2m v qV 



(A 3 (q 2 )-A (q 2 )) 



(25) 



In the rest-frame of the heavy meson P, the product of 
the four-momentum of the meson P and the polarization 
vector of the photon satifies pp-e = 0. Then the last term 
in the above equation can be dropped. With py = p~ n 



5 



we obtaine 
Ald = 



where: 



eG F V Qq Q q 



v 



V2 

■ LD^ap/icrP 



uiY 1 (1 - 75) Vp 



x {Wld 



V LD 



A v - 



C q vfv 



-m,y + irnv^v mp + my 
e^Cqyfv mp + m v 
-my + im v T v pp-p 1 



V(q 2 ) 



IV Numerical Result and Discussion 



(26) 



(27) 
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The numerical calculation is performed in the center- 
of-mass frame of the heavy meson, and the momentum 
of the photon is taken as (.Ely, 0, 0, — -Ely). For the input 
parameters, the masses of the quarks are taken as 

m u = m d = 0.08GeV, m s = 0.30GeV, 
m b = 4.98GeV, m c = 1.54GeV, 

which are taken to be consistent with that used to derive 
the wave function of the heavy meson in the relativistic 
potential model in Ref. |6|. 

The electron and neutrinos are taken to be massless, 
the masses of the other leptons are taken from PDG 
0. The Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements are: 

V cd = 0.2259, V cs = 0.974, V ub = 0.00389. 

The infrared divergence appears as the energy of the 
real photon goes to soft limit or the momentum of the 
photon is parallel to the momentum of a massless lepton. 
This divergence can be canceled when the decay width 
of the radiative leptonic decay is added with the pure 
lcptonic decay width, in which one-loop radiative correc- 
tions are included [17[ . The radiative leptonic decay with 
the energy of the photon lower than the experimental res- 
olution can not be distinguished from the pure leptonic 
decay. Only photons with the energy larger than the 
experimental resolution can be distinguished. Therefore 
the radiative leptonic decay width depends on the pho- 
ton energy resolution. The photon energy resolution can 
be a few MeV in experiment [l8| . The dependence of the 
decay width on the resolution AE 7 is shown in Table U 
and Fig. [3] For example, if taking AE~ t = 10 MeV, the 
decay width and branching ratio of D — > "feD e are 

T(D je9 e ) = 1.98 x 10~ 18 GeV, 

« (28) 
Br(D -> 7 eP e ) = 3.29 x 10~ 6 . 

In the following all the numerical calculation is performed 
by taking the resolution Ai? 7 = 10 MeV. 



FIG. 3: Decay widths of D — > 'yev as a function of AE- 



TABLE I: The dependence of the decay width of D —± jeu 
on the photon resolution AE-,. 



AEj T D ~^ lW AE-, r D -_ >7ii; 
(MeV) (xlCT 18 GeV) (MeV) (xlO~ 18 GeV) 



5 


2.30 


55 


1.25 


10 


1.98 


60 


1.22 


15 


1.80 


65 


1.19 


20 


1.67 


70 


1.16 


25 


1.57 


75 


1.14 


30 


1.49 


80 


1.12 


35 


1.43 


85 


1.10 


40 


1.38 


90 


1.08 


45 


1.33 


95 


1.06 


50 


1.29 


100 


1.04 



We also checked that the decay width is not sensitive 
to the cutoff of the angle between the momentum of the 
photon and the electron. We find taking different val- 
ues for the angle cutoff A9, the result changes extremely 
slowly. 

To show the contribution of each diagram in Fig[TJ we 
list each diagram's contribution to the decay width in 
Table HU It is shown that, although the amplitude of the 
diagram with the photon emitted from the heavy quark 
is suppressed by a factor of the heavy quark mass in the 
dominator of the propagator, which is like . _J _ m , for 

the cases of D and D s mesons, the mass of c quark is not 
large enough, so the contributions of diagrams a, b and 
c are all at the same order, but for the case of B decay, 
the suppression is large, the contribution of diagram b 
dominates. It can also be shown that, the contributions 
of the diagrams in Fig JT] interfere destructively, especially 
in the case of D and D s mesons, this is consistent with 
Ref. [!]. 

We present the branching ratios with only the short- 
distance contributions for all the decay modes in TA- 



6 



TABLE II: The decay width contributed by each diagram in 
Fig. \T\ In the first column, r a , Ft and r c means the contri- 
butions of diagrams a, b and c, respectively, T a +b the contri- 
bution of the diagrams a and b, T a+ t+ c the total contribution 
of diagrams a, b and c. 



Width D~ 
(GeV) 



D- 



B~ 



r a 


1.37 x 10" 


17 


3.06 x 10" 


16 


9.14 X 10" 


21 


r b 


1.35 x 10" 


17 


1.94 x 10" 


16 


1.01 X 10" 


18 


r c 


9.78 x 10" 


18 


3.23 x 10" 


16 


5.83 x 10" 


20 




9.86 x 10" 


18 


3.08 x 10" 


16 


0.96 x 10" 


18 




1.98 x 10" 


18 


8.46 x 10" 


18 


0.82 x 10" 


18 



TABLE III: Branching ratios of radiative leptonic decays of 
B and D mesons. BRsd is the branching ratios with only 
the short-distance contribution. BRld stands for branch- 
ing ratio with only long-distance contribution. For D and B 
mesons, BRldi is the branching ratios of long-distance con- 
tribution via p meson, while for D s meson, the long-distance 
contribution is via <f> meson. BRld2 is the branching ratios 
of long-distance contribution via cj meson. 



Modes 


BRsd 




BRldi 




BRlD2 




D~ - 


-> 7e^ e 


3.3 x 10" 


(i 


7.5 x 10" 


(i 


6.3 x 10" 


6 


D~ - 


-» 7M^m 


1.6 x 10" 


5 


7.3 x 10" 


6 


6.1 x 10" 


6 


D' - 


-¥ ^tu t 


5.5 x 10" 


9 


9.1 x 10" 


10 


7.6 x 10" 


10 


Ds - 


->• -yeue 


6.8 x 10" 


(i 


1.0 x 10" 


4 






D~ - 


-> 7M^m 


2.0 x 10" 


4 


1.0 x 10" 


4 






D7 - 


-¥ ^tv t 


1.1 x 10" 


6 


6.5 x 10" 


8 






B~ - 


-» 7£^ e 


2.1 x 10" 


(i 


5.8 x 10" 


7 


4.0 x 10" 


7 


B~ - 


+ 7/^ 


2.0 x 10" 


6 


5.8 x 10" 


7 


4.0 x 10" 


7 


B~ - 




1.6 x 10" 


6 


3.8 x 10" 


7 


2.7 x 10" 


7 



BLE IIIII The branching ratio for D~ —> jtD t is very 
small, because the mass of r is very large, the phase 
space for this decay mode is too small. 

The short-distance branching ratios obtained in this 
work is slightly smaller than the previous works Q • More 
details about the quark momentum distribution are in- 
cluded in this work by using wave function obtained in 
the relativistic potential model. 

To calculate the long-distance contribution, the tran- 
sition amplitude V — > 7 is needed. The transition ampli- 
tude is related to the decay constant defined in eq. (l2l"1) . 
Using the data on the decay rate of V —> l + l~ given in 
PDG 0, the decay constants of the vector mesons can 
be extracted as 



f p = 0.169869GeV, 
f u = 0.154631GeV, 
U = 0.231784GeV. 



(29) 



The q 2 dependence of the form factors defined in the 
hadronic matrix element < V\ Jy_ A \P > are taken to be 
the usual pole approximation as 



V(q 2 



-, A^q 2 



The parameters in the form factors for D — > V and B — > 
V transitions are quoted from Refs. [l9| and (20j . they 
are 



h v D ^ p = 1-0, 

M Vo ^ p = 2.5 GeV; 
hv B ^ p = 0.323, 
h a i B ^ p = 0.242, 
h VB ^ = 0.293, 
h a i B ^ v = 0.219, 



M AlL 



= 0.5; 

p = 2.5 GeV; 
= 38.34 GeV 2 ; 
p = 37.51 GeV 2 ; 
"= 37.45 GeV 2 ; 
= 37.01 GeV 2 . 



For the form factor of D — > w, we assume it is the same 
as that of D — s- p. The form factors of D s — > <fi is taken 
from 211 as 



1.21, h„ 



0.55, M- 



2.08 GeV. 



Mi 



The form factor A\ for D s — > </> transition is approxi- 
mated as a constant because the q 2 dependence of Al is 
very weak 21 1 . 

With the parameters given above, the long-distance 
contributions to the decay width for each decay mode 
can be estimated, they are listed in Table IIII| where 
the short-distance and long-distance contributions to the 
branching ratios of each decay mode are presented sep- 
arately . It shows that for decays of D and D s mesons, 
long-distance contributions are as large as or even larger 
than the short-distance contributions, while for the case 
of B decays, short-distance contributions dominate, long- 
distance contributions are roughly 4 ~ 5 times smaller 
than short-distance contributions. 

To get the total decay width, including both the short 
and long-distance contributions, one has to know the 
relative phase between the long and short-distance am- 
plitudes. Unfortunately we do not konw the relative 
phases exactly upto now. We have to leave the rela- 
tive phases as free parameters. To show how the decay 
width depends on the relative phases, we show the de- 
cay widths of B, D and D s — > -feD e decays in Fig. 0] 
as an example. In the case of D — > ^ev e decay, be- 
cause the long-distance contributions are as important 
as the short-distance contributions, the relative phases 
between the long and short-diatance contributions can 
affect the decay widths considerably, the decay widths 
can change several times. For D s — > 7ez/ e decay, the 
long-distance contribution dominates (see Table Mil) , the 
dependence of the total branching ratio on the relative 
phase is weak. While for B — > jeD e meson decay, the 
amplitudes of the short and long-distance contributions 
are at the same order, therefore the decay width still de- 
pends on the relative phase severely. For the case of the 
other decay modes, the dependence of the total branch- 
ing ratios on the relative phases are illustrated numeri- 
cally in Table I IV! The situation is similar to the B, D 
and D s — > ^ev e decay modes. The contribution of the 
long-distance physics is important, in general the branch- 
ing ratio heavily depends on the relative phase between 
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-2 -] I 2 3 



FIG. 4: From up to down, they are decay widths (GeV) of 
D~, B~ and D~ — > ^e~v e decays, including both the short 
and long-distance contributions, as functions of <j) p and tfru, in 
the cases of D~ and B~ meson decays, and 4><i> m the case of 
meson decay. 



the long and short-distance contributions. Some decay 
modes can be greatly enhanced by the long-distance con- 
tributions. The branching ratio of D s — » ^ev e decay can 
be enhanced from the order 10~ 6 to 10~ 4 . 

V Summary 

The radiative leptonic decays of B, D and D s — > 
are studied in this work. The short-distance contribution 
is calculated by using the wave functions of the heavy 
mesons obtained in the relativistic potential model, more 
details about the quark-momentum distribution are in- 



cluded in this work. In addition to the short-distance 
contribution, the long-distance contribution is also esti- 
mated in the vector meson dominance model. The study 
shows that the long-distance contributions can seriously 
affect the decay rates. The branching ratio of D~ — > ^ev e 
can be enhanced to the order of 10 -4 , which should only 
be at the order of 10~ 6 if only considering the short- 
distance contribution. 
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TABLE IV: Total branching ratios of the radiative leptonic decays including both the short and long-distance contributions. 
In case of D meson decays, for illustration, the relative phase cf> = <j> w and cf> p — are taken. While for B meson decays, 
4> = 4> p = 4> w is taken, and in case of D s meson decays, is cj)^. 



Modes 


BR t ot(4> = 


0"' 


) BR tot ((j) = 


30°) BRtot(4> = 


60°) BR tot (ct> = 


90°) 


BR t ot(<t> = 120°) BR t ot(4> = 150°) 


D — > 'yeve 


4.U X 1U 


■b 


1.1 X 1U 




3.0 x 10 


— b 


4.1 x 10 




3.4 x 10" b 


1.5 x 10" b 


D -¥ 7^ii? M 


2.5 x 10" 


■5 


3.1 x 10 


— 5 


4.8 x 10 


-5 


5.9 x 10 


-5 


5.2 x 10~ 5 


3.5 x 10" 5 


D — > "1TU T 


5.7 x 10 - 


■9 


7.5 x 10" 


-9 


1.2 x 10" 


-8 


1.4 x 10" 


-8 


1.2 x 10~ 8 


8.0 x 10" 9 


B — y ^eve 


6.7 x 10" 


■6 


8.9 x 10 


-b 


1.3 x 10" 


-b 


1.4 x 10" 


-b 


1.2 x 10"° 


8.2 x 10"° 


B -»■ 7/i!? M 


6.8 x 10" 


6 


8.9 x 10 


-6 


1.3 x 10 


-5 


1.4 x 10 


-5 


1.2 x 10~ 5 


8.3 x 10 -6 


B — > ^TV T 


5.3 x 10" 


6 


6.7 x 10 


-6 


9.2 x 10" 


-6 


1.0 x 10" 


-5 


8.8 x 10~ 6 


6.3 x 10" (i 


D 3 — > 7ei/ e 


1.5 x 10" 


■4 


1.3 x 10 


-4 


0.9 x 10" 


-4 


0.7 x 10" 


-4 


0.9 x 10" 4 


1.3 x 10" 4 


D s ->■ 7^ 


4.5 x 10" 


■4 


4.4 x 10 


-4 


4.1 x 10 


-4 


3.9 x 10 


-4 


4.1 x 10~ 4 


4.4 x 10" 4 


D s — > -yru T 


1.0 x 10" 


6 


1.0 x 10 


-6 


1.2 x 10 


-6 


1.3 x 10 


-6 


1.2 x 10 -6 


1.0 x 10" 6 



